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Abstract 

The dual-tree complex wavelet transform (DT-CWT) is known to ex- 
hibit better sliift-invariancc than the conventional discrete wavelet trans- 
form. We propose an amplitude-phase representation of the DT-CWT 
which, among other things, offers a direct explanation for the improve- 
ment in the shift-invariancc. The representation is based on the shifting 
action of the group of fractional Hilbert transform (fHT) operators, which 
extends the notion of arbitrary phase-shifts from sinusoids to finite-energy 
signals (wavelets in particular). In particular, we characterize the shifta- 
bility of the DT-CWT in terms of the shifting property of the fHTs. At 
the heart of the representation are certain fundamental invariances of the 
fHT group, namely that of translation, dilation, and norm, which play a 
decisive role in establishing the key properties of the transform. It turns 
out that these fundamental invariances are exclusive to this group. 

Next, by introducing a generalization of the Bedrosian theorem for the 
fHT operator, we derive an explicitly understanding of the shifting action 
of the fHT for the particular family of wavelets obtained through the mod- 
ulation of lowpass functions (e.g., the Shannon and Gabor wavelet). This, 
in effect, links the corresponding dual-tree transform with the framework 
of windowed- Fourier analysis. 

Finally, we extend these ideas to the multi-dimensional setting by in- 
troducing a directional extension of the fHT, the fractional directional 
Hilbert transform. In particular, we derive a signal representation in- 
volving the superposition of direction-selective wavelets with appropri- 
ate phase-shifts, which helps explain the improved shift-invariance of the 
transform along certain preferential directions. 

1 INTRODUCTION 

The dual-tree complex wavelet transform (DT-CWT) is an enhancement of 
the conventional discrete wavelet transform (DWT) that has gained increas- 
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ing popularity as a signal processing tool. The transform, originally proposed 
by Kingsbury [12] to circumvent the shift-variance problem of the decimated 
DWT, involves two parallel DWT channels with the corresponding wavelets 
forming approximate Hilbcrt transform pairs [16]. Wc refer the reader to the 
excellent tutorial [16] on the design and application of the DT-CWT. 

In this contribution, we characterize the dual-tree transform from a com- 
plementary perspective by formally linking the multircsolution framework of 
wavelets with the amplitude-phase representation of Fourier analysis. The lat- 
ter provides an efficient way of encoding the relative location of information in 
signals through the phase function that has a straightforward interpretation. 
Specifically, consider the Fourier expansion of a finite-energy signal f{x) on 
[0,L]: 

f{x) — ao + ai cos{ujox) + a2 cos(2ujox) + ■ ■ ■ 

+ bi sin(woa;) -|- 62 sin(2a;oa^) (1) 

Here loq {ujqL = 27r) denotes the fundamental frequency, and ao, ai, 02, • • ■■, and 
61,62, .. . are the (real) Fourier coefficients corresponding to the even and odd 

harmonicsrespectively. Now, by introducing the complex Fourier coefficients 
Cn = o^n ~jbn ^nd by expressing them in the polar form c„ = |c„|e-''^", ^ 0„ < 
27r, one can rewrite (1) as 

00 00 
f{x) = I c„ I f cos (j)n cos{nu)ox) - sin0„sin(na;oa;)) = ^ |c„|(p„(a: -h -r„) 

n=0 n=0 

(2) 

with Tn = (pn/nuiQ specifying the displacement of the reference sinusoid (Pn{x) = 
cos{nLU()x) relative to its fundamental period [0,L/n\. The above amplitude- 
phase representation highlights a fundamental attribute of the shift parameter 
f„: it corresponds to the shift f that maximizes |(/(-)) + 'f'))\, the correla- 
tion of the signal with the reference cp„(a;). The corresponding amplitude |c„| 
measures the strength of the correlation. 

As far as signals with isolated singularities (e.g. piecewise-smooth signals) 
arc concerned, the wavelet representation employing the dilated-translated 
copies of a fast-decaying oscillating waveform - has proven to be more efficient. 
Moreover, the added aspect of multiscale representation allows one to zoom onto 
signal features at different spatial resolutions. Complex wavelets, derived via 
the combination of non-redundant wavelet bases, provide an attractive means 
of recovering the crucial phase information. In particular, the phase relation 
between the components (of the complex wavelet) is used to encode the relative 
signal displacement (besides offering robustness to interference). The DT-CWT 
is a particular instance where the components are related via the Hilbert trans- 
form [16, 15]. 
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1.1 Main Results 



Analogous to the fact that the complex Fourier coefficients in (2) are derived 
from the (primitive) analytic signals e^"""^ = Vn{x) + jJif(Pn{x): 



the DT-CWT coefficients are obtained by projection the signal on to the dilated- 
translated copies of the analytic wavelet \E'(a;) = tpi^) + j-^i^ix). The cen- 
tral idea of this paper is the identification of the wavelet counterparts of the 
phase-shifted sinusoids (p„(a; + T„) in (2). In particular, the died shifted 

wavelets are derived by the action of the single-parameter family of fHT oper- 
ators {^IreR (for definition and properties see §2) on the reference wavelet 
■^(2;). The shift parameter r controls the shifting action il){x) ^ :f{rtp{x) of the 
fHT, and, in effect, resulting in the realization of a continuously-defined fam- 
ily of shifted wavelets. This action has an important connotation in relation 
to pure sinusoids: ^(p„(x) = cos(na;oa; + ttt). In fact, the amplitude-phase 
representation of the DT-CWT is derived in §3.1 by generalizing the following 
equivalent expression of (2): 



which, in turn, is based on the above-mentioned phase-shift action of the fHTs. 
The significance of either representation is that they allow us to provide a precise 
characterization of the shiftability of the associated reference functions in terms 
of the shifting action of the fHT. 

Motivated by this connection, we make a detailed study of the group of fHT 
operators in §2.1. In particular, we highlight their invariance to translations and 
dilations which allows us to seamlessly integrate them into the multiresolution 
framework of wavelets. Moreover, the observation that the above-mentioned 
invariances are exclusively enjoyed by the fHT group (cf. Theorem 3.1) makes 
the shiftability of the dual-tree transform unique. For the particular family of 
dual-tree transforms involving HT pairs of modulated wavelets, we derive an 
explicit characterization of the shifting action of the fHT in §3.2. If the dual- 
tree wavelet is not modulated, we can still characterize the action of the fHT by 
studying the family of fractionally-shifted wavelets {^V'ItsR) and we do this 
explicitly for the particular case of spline wavelets in §5. Finally, we extend the 
proposed representation to the multi-dimensional setting in §4 by introducing 
certain directional extensions of the fHT. 

The above results have certain practical iniplic;ations. In §3.3, we propose 
certain measures for accessing the quality of the factors, namely the HT cor- 
respondence and modulation criterion, that are fundamental to the shiftability 
property of the dual-tree wavelets. These metrics could prove useful in the 
design of dual-tree wavelets with a good shift-invariance property. 



Cn = an-jbn = {f{x),e- 






n=0 
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2 THE FRACTIONAL HIBERT TRANSFORM 



In what follows, the Fourier transform of a function f{x) defined over R"^ {d > 1) 
is specified by f{uj) = J-^^ exp {—jiJ^ x)dx, where uP^x denotes the usual 
inner-product on R''. The other transform that plays a significant role is the 
Hilbert transform (HT) [6, 18] ; we will denote it by . In particular, we shall 
frequently invoke the Fourier equivalence 

Jf'/(x)^-isign(a;)/(a;) (3) 

characterizing the action of the HT on L^(R), the c;lass of finite-energy signals"^. 
Three fundamental properties of the HT that follow from (3) are its invariance to 
translations and dilations, and its unitary (norm-preserving) nature. Moreover, 
we shall use to denote the identity operator {,f f){x) = f{x). 

We begin with a detailed exposition of the relevant characteristics of the fHT 
that forms the cornerstone of the subsequent discussion. There exit several def- 
initions of the fHT in the signal processing and optics literature [23, 19, 13, 9]; 
however, for reasons that will be obvious in the sequel, we propose to formulate 
it as an interpolation of the "quadrature" identity and HT operator using con- 
jugate trigonometric; func;tions. In particular, we define the fHT operator H^, 
corresponding to the real-valued shift parameter r, as 

= cos(7rT) .y - sin(7rT) J^. (4) 

This definition is equivalent to the formulation introduced in [13, 9], but 
differs from the ones in [23, 19] up to a complex chirp. The important aspect 
of the above operator-based formulation is that it directly relates the fHT and 
its properties to the more fundamental identity and HT operator, which are 
identified a posteriori as special instances of the fHT: ^ = and M' = ^-xji- 
In view of (4), we would like to make note of fact that ^ and H have a non-linear 
correspondence: 

Jf 2 = -J, (5) 
and act in "quadrature" in the sense that 

W,Jf/)=0 (/eL2(R)). (6) 

These come as a direct consequence of definition (3) and certain properties of 

the inner- product. 

As far as the domain of definition of (4) is concerned, note that both J' and 
^ act as bounded operators (with a bounded inverse) on L^'(R) for 1 < p < oo 
[4], and so does ^Ht. In particular, the fHT admits the following equivalent 
specification on L-^(R): 

Hrf{x) ^ exp (jTTT sign(u;)) /(u;) (7) 

that comes as a consequence of equivalence (3). We shall henceforth invoke (4) 
and (7) interchangeably in the context of finite-energy signals. 

^ The domain can also be extended to include distributions such as the Dirac delta and the 
sinusoid [4, Chapter 2] 



4 



2.1 Chciracterization of the fHT 



As remarked earlier, most of the characteristic features of the constituent iden- 
tity and HT operators carry over to the family of fHT operators {^jreR- 
particular, the following properties of the fHT can be readily derived: 
(PI) Translation- invariance: {Hrf){x — y) = ^{/(- — y)}{x) for all real y. 
(P2) Dilation-invariance: (5^/)(Ax) = :?^-{/(A-)}(a;) for all positive A. 
(P3) Unitary nature: {y{rf,!f{rg) = {f,g); in particular, = ||/|| for all 

/,gin L2(R). 

(P4) Composition law: y{ri^2 = ^i+t2- 

(P5) Phase-shift operator: :Kr{cos{uiox)} = cos(a;oa^ + ttt). 

Indeed, (PI) and (P2) are immediate consequences of the dilation- and 
translation-invariance of and Jif; (P3) follows from the unitary frequency 
response (7) and Parseval's identity; and (P4) follows from (5): 

i^^i^^ = (cos(7rTi) y - sin(7rTi) Jf ) ( cos(7rT2) - sin(7rr2) Jf) 



Finally, it is the quadrature-shift action cos(a;oa;) i— > sin(a;oa;) of Jif that results 



^ cos(a;oa;) = cos(7rr) cos(a;o2;) — sin(7rr) J^{cos(a;o2;)} = cos(a;oa; -|- ttt); (8) 

this is also justified by the frequency response of the fHT. 

As will be discussed shortly, properties (PI), (P2) and (P3) play a crucial 
role in connection with wavelets. The composition law (P4) tells us that the 
family of fHT operators is closed with respect to composition. Moreover, as the 
identity Hr^-r = !Ho = ^ suggests, the inverse^ fHT operator is also a fHT 
operator specified by = These closure properties can be summarized 
by following geometric characterization: 

Proposition 2.1 The family of fHT operators forms a commutative group on 
LP(R),l<p<oo. 

In this respect, note the marked resemblance between the family of fHT opera- 
tors and the commutative group of translation operators that play a fundamen- 
tal role in Fourier analysis; the relevance of the former group in connection with 
the DT-CWT will be demonstrated in the sequel. Moreover, the finite subgroup 
{j^,^,—J^,—y{} of self-adjoint"^ operators (see Fig. 1) is worth identifying. It 
is the smallest subgroup containing the in-phase/quadrature operators that play 
a fundamental role in the dual-tree transform. 

■^Note that, for a given t, there exists infinitely many t' such that identity 9{r^i = 
holds; this comes as a consequence of the periodicity of the trigonometric functions involved 
in the definition. One can easily factor out the periodic structure by identifying ri and T2 iff 
Hr-^ = ; in effect, this equivalence relation results in the specification of equivalence classes 
of fHTs. However, for the simplicity of notation, we shall henceforth use Ht to denote both 
the equivalence class and its representatives. 

^self-adjoint up to a sign: T* = ±T for each T in the subgroup. 



= COS 




in (P5): 
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Figure 1: Geometrical interpretation of the continuous fHT group using the 
isomorphic unit-circle group (the multiplicative group of complex numbers 
having unit modulus); the correspondence is < (cos(7rr), sin(7rT)). 

2.2 The Wavelet Context 

Similar to the Hilbert transform, the fHT perfectly fits the wavelet framework. 
The implication of properties (PI) and (P2) is that the fHT of simultaneous 
dilates and translates of a wavelet is a wavelet, dilated and translated by the 
same amount. This has fundamental ramifications in connection with dyadic 
wavelet bases generated via the dilations and translations of a single mother- 
wavelet ipi^x). In particular, let 4'i,k{x) and Si^k denote the dilated-translated 
wavelets Si_fcV(2;) = 2'/^i/)(2*x — k), {i, k) G 7? and the corresponding (normal- 
ized) dilation-translation operators, respectively. Then the commutativity 

Hr Ei^k = Ei,k^r (9) 

holds for all real r and integers i and k. The significance of (9) is that it allows 
us to conveniently factor out the pervasive dilation-translation structure while 
analyzing the action of fHTs on wavelet bases. 

On the other hand, a fundamental consequence of the isometry property 
(P3) is that maps a Riesz basis onto a Riesz basis; in particular, if {ijji^k} 
forms a wavelet basis of L^(R), then so does {y^rfpi.k}- In fact, J^r preserves 
biorthogonality: if {"01, fc} and {V'i'.fe'} constitute a biorthogonal wavelet basis 
satisfying the duality criteria {ipi,k,4^i',k') = S[i — i']S[k — k'], then also we have 
that 

{'HripiM.'Hr^i'M') = {ipi,k,'ipi'.k') = S[i - i']S[k - fc'], 
signifying that {Hripi^k} and {7irV'i',fe'} form a biorthogonal basis as well. 



6 



3 SHIFTABILITY OF THE DUAL-TREE TRANS- 
FORM 



3.1 Multiscale Amplitude-Phase Representation 

We now derive the amplitude-phase representation of the DT-CWT based on 

the shifting action of the fHT. As remarked earlier, the parallel is grounded 
on the observation that instead of the quadrature sinusoids {cos(na;oa;)} and 
{sin(na;oa^)}) the DT-CWT employs two parallel wavelet bases, {tpi,k} and 
{''Pi k}' derived via the dilations-translations of the wavelets ijj{x) and %i''{x) that 
form a HT pair: tp'i^) = J^i^ix). A signal f{x) in L^(R) is then simultaneously 
analyzed in terms of these wavelet bases yielding the wavelet expansions 

The analysis coefficients [k] and bi [k] are specified by the projections onto the 
dual wavelet bases and {V'^ j,}: 

= (/, V^i.fe), and bi[k] = {f,ip-j,). 

The dual wavelet bases are implicitly related to the corresponding primal bases 
through the duality criteria (V'i,fc, V'p.n) = S[i — p,k — n] and (V'- fc, V'p,™) = 
6[i — p, k — n]. A fundamental consequence of the unitary property of the HT is 
that these dual bases can be generated through the dilations-translations of two 
dual wavelets, say tl>{x) and tl>'{x), that form a HT pair as well: V''(a;) = M''ij){x) 
[7] . In particular, by introducing the complex wavelet 

^{x) = \{'^{x)+j^'{x)) 

and its dilated-translated versions ^i^k{x) - the analytic counterpart of the 
complex sinusoids exp(jna;ox) - the dual-tree analysis can simply be viewed 
as the sequence of transformations f{x) i— > Ci[k] = {f,^i^k) resulting in the 
complex analysis coefficients {ci[k]}(^i^k)ez'^- 

Our objective is to derive a representation of f{x) in terms of the modulus- 
phase information Ci[k] = |ci[fc]|e^'^'['^', and the reference wavelets tpij.{x). In 
particular, by combining the expansions in (10) and by invoking (9), we arrive 
at the following representation: 

E {ai[k]iJi,k{x) + bi[k]^lk{x)) 

(i,fe)6Z2 

= E %[/c]A{V'*,/c(a;)} 
(i,fe)ez2 

= h[k]\^i,k{^{x;Ti[k])}, (11) 

(i,fe)ez2 
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where the synthesis wavelet tp(x; Ti[k]) = H^^^^ip{x) is derived from the mother 
wavelet ^(x) through the action of the fHT corresponding to the shift Ti[k] = 
(j>i[}t\/-K. The above multiresolution amplitude-phase representation provides 
two important insights into the signal transformation 

The first of these is derived from the observation that the fractionally-shifted 
wavelets tp{x]Ti[k]) in (11) play a role analogous to the phase-shifted sinusoids 
(p„(a; + T„) in (2). In particular, (11) offers a rigorous interpretation of the 
ampitude-phase factors: while \ci[k]\ indicates the strength of wavelet correla- 
tion, the relative signal displacement gets encoded in the shift Ti[k] correspond- 
ing to the most "appropriate" wavelet within the family {i^'^i'i,fe}TeR- 

The "shiftable-wavelet representation" in (11) also offers an explanation for 
the improved shift-invariance of the dual-tree transform which is complementary 
to the frequency-domain argument given by Kingsbury in [12]. It is well-known 
that the signal representation associated with the (critically-sampled) discrete 
wavelet transform is not shift-invariant; in particular, the uniform sampling 
of the translation-parameter of the continuous transform limits the degree of 
shift-invariance of the transform. The fact that the over-sampled dual-tree 
transform tends to exhibit better shift-invariance can then be explained in terms 
of the associated shiftability of the transform. Indeed, the fractional-shifts of 
the reference wavelets around their discrete translates partially compensates for 
the limited freedom of translation. 

It is clear that the invariances of the fHT group were central to the derivation 
of (11). The following result establishes the fHTs as the only complete family 
of operators that exhibits such characteristics: 

Theorem 3.1 (Uniqueness of the fHT) 

A unitary linear operator T on (R) is invariant to translations and dila- 
tions if and only if it can he represented as 

T = cose J' -sine (12) 

for some unique 6 in (— 7r,7r]. 

The above result (proof provided in Appendix 7.1) signifies that any family of 
unitary operators that simultaneously commutes with translations and dilations 
is isomorphic to the fHT group. This provides significant insight into the rep- 
resentation in (11) since it is these fundamental invariances that facilitate the 
incorporation of the fHT into the wavelet framework. 

As discussed next, it turns out that the shifted wavelets can be very explicitly 
characterized for certain classes of wavelets which provides a deeper insight into 
the above signal representation. 
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Figure 2: Idealized spectrum of a modulated wavelet. The spectrum has pass- 
bands over uq — ^ < \u)\ < loq + with local axes of symmetry at a; = ±loo- 

3.2 Modulated Wavelets: Windowed-Fourier-like Repre- 
sentation 

A wavelet is, by construction, a bandpass function. Specifically, if the wavelet 
is a modulated function of the form 

ip{x) = (p{x) cos (ujox + ^o) , (13) 

where (p{x) is bandlimited to [—17,17] (for some arbitrary 17), and > 17 [cf. 
Fig. 2], we can make precise statements on the representation (11) correspond- 
ing to the dual-tr(K^ transform involving such a modulated wavelet and its HT 
pair. In order to do so, we provide a very generic result that allows us to extend 
the phase-action in (8) to modulated functions: 

Theorem 3.2 (Generalized Bedrosian Identity) 

Let f{x) and g{x) he t,wo real-valued functions such that the support of /(lu) 
is restricted to (—17,17), and that g{uj) vanishes for |a>| < 17 for some arbitrary 
frequency 17. Then the fHT of high-pass function completely determines the fHT 
of the product: 

H,{f{x)g{x)} = f{x)^g{x). (14) 

Informally, the above result (cf. Appendix §7.2 for a proof) asserts that 
the fHT of the product of a lowpass signal and a highpass signal (with non- 
overlapping spectra) factors into the product of the lowpass signal and the fHT 
of the highpass signal. Note that, as a particular instance of Theorem 3.2 
corresponding to r = —1/2, we recover the result of Bedrosian [3] for the HT 
operator. An important consequence of Theorem (3.2) is that, for a wavelet is 
the form (p{x) cos{u)ox) with (p{x) is bandlimited to {—uo,u)o), the fHT acts on 
the phase of the modulating sinusoid while preserving the lowpass envelope: 

^{(p{x)cos{u)ox)} = (p{x)cos{coox + ttt). (15) 
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In particular, the above modulation action allows us to rewrite the signal rep- 
resentation in (11) as 

fixed window variable amp— piiase oscillation 

^ (Pi,fc(x) Ei^k[ \ci[k]\ COS {ujox + ^o + 'KTi[k])Y (16) 
(i,fc)ez2 

where ^)i^k{x) = 'E.i k(p(x) denotes the (fixed) window at scale i and translation 
k. This provides an explicit interpretation of the parameter [k] as the phase- 
shift applied to the modulating sinusoid of the wavelet. In this regard, its role is 
therefore similar to that of the shift parameter in in the Fourier representation 
(2). In effect, while the localization window (pi.k{x) is kept fixed, the oscillation 
is shifted to best-fit the underlying signal singularities/transitions. In this light, 
one can interpret the associated dual-tree analysis as a multiresolution form of 
the windowed-Fourier analysis, with the fundamental ciifferencx; that, instead of 
analyzing the signal at different frequencies, it resolves the signal over difi^erent 
scales (or resolutions). 

Two concrete instances of such modulated wavelets are: 

• Shannon wavelet: The Shannon wavelet is constructed from the Shan- 
non multiresolution [14]; it is specified as 

*W = si„c(^)»«(?^i^). 

and it dilates-translates constitute an orthonormal wavelet basis of L^R. 
Many wavelet families converge to the Shannon multiresolution as the 
order increases [1]; e.g., the orthonormal Battle-Lemarie wavelets [2], and 
the interpolating Dubuc-Deslauriers wavelets [10]. 

• Gabor wavelet: The sinc(.x) envelope of the Shannon wavelet results 
in an "ideal" frequency resolution but only at the expense of poor spa- 
tial decay. As against this, wavelets modeled on the Gabor functions [11] 
(modulated Gaussians) exhibit better space-frequency localization. More- 
over, as in the case of the Shannon wavelet, they are quite generic in 
nature since several wavelet families closely resemble the Gabor function. 
For example, the B-splines wavelets (a semi-orthogonal family of spline 
wavelets) asymptotically converge to the Gabor wavelet 

4!{x) = g"{x) cos {ojox + ^o), 

where g°'{x) is a Gaussian window that is completely determined by the 
degree of the spline [20]. In fact, based on this observation, a multiresolu- 
tion Gabor-like transform was realized in [7] within the framework of the 
dual-tree transform. This Gabor-like transform involved the computation 
of the projections 

fix) ^ {f{x),Ei,k^{x)) {i, fc e Z) (17) 
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on to the dilatcs-translatcs of the Gabor-hkc wavelet ^{x) = tjj{x) + 
jJifijj{x), and was realized using the usual dual-tree transform correspond- 
ing to the spline wavelets 4>{x) and J^^p{x). We would, however, like to 
point out that the representation in (16) corresponds to a situation where 
the role of the analysis and synthesis functions have been reversed, namely 
one in which the signal is analyzed using the dual complex wavelet ^{x), 
and where the Gabor-like wavelet ^'(a;) is used for reconstruction. 

Graphical Illustrations: Fig. 3 shows quadrature pairs (:^V(2^)i^-i-i/2V'(^)) 
of Shannon- like (resp. Gabor-like) wavelets corresponding to different r. The 
fHT pair {^,^+1/2) drives the modulating oscillation to a relative quadra- 
ture that are localized within a common sinc-like (resp. Gaussian-like) window 
specified by \^'tp{x) + j^^i/2tp{x)\. 

In figure 4, we demonstrate the shiftability of the dual-tree transform using 
a Gabor-like B-spline wavelet ip{x) of degree 3 as the reference. To this end, 
we consider the step input f{x) = sign(a; — xq) that has a discontinuity at xq. 
The M-level decomposition of this signal (cf. [7] for implementation details) in 
terms of the conventional DWT is given by 

M 

e i=l k 

with (PM,e{x) denoting the translates of the coarse representation of the scaling 
function. As for the DT-CWT, we have the representation 

M 

fi^) = '^PiVM,e{x) + '^PWM,n{x) + '^'^\ci[k]\ il)i,k{x;Ti[k]). (19) 

i n i—1 k 

The idea here is to demonstrate that the shifted wavelets in (19) respond better 
to the signal transition than the wavelets in (18); that is, the oscillations of the 
shifted wavelets have a better lock on the singularity at Xq. Figure 4 shows the 
reference wavelet ^i,k{x) and the shifted wavelet fpi,k{x;Ti[k]) corresponding to 
a specific resolution i = J and translation k = [xo/'2'^] (around the position 
of the singularity at the coarser resolution) . Also shown in the figures arc the 
step input and the fixed Gaussian-like localization window of the wavelet. The 
oscillation of the shifted wavelet is clearly seen to have a better lock on the 
transition than the reference. The magnitude of the signal correlation in either 
case justifies this observation as well. 

3.3 Quality Metrics for Dual- Tree Wavelets 

The shiftability of the DT-CWT was established based on two fundamental 

properties of the wavelets, ipi{x) and ijj2{x), of the two branches: 
(CI) HT correspondence: ip2{x) = Jf^pi{x), 

(C2) Their modulated forms: tpi{x) = (p{x) cos(woa;-|-^o)) and il>2{x) = (p{x) sin(a;oa;-|- 
^0), with the support of (p{oj) restricted to [—uo,iOo]. 
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(1a)T= -1 



(1b) T = 



(1c) T = 1/2 




Figure 3: Quadrature pairs of orthonormal spline (resp. B-spline) wavelets re- 
sembling the Shannon (resp. Gabor) wavelet: Solid (blue) graph: :Hrip^{x); 
Broken (red) graph: yir+i/2''P^i^)': ^^nd Solid (black) graph: Common localiza- 
tion window given by \^'iIj^{x) + jy-Cr+i/2''P^{x)\- 



Table 1: Quality metrics for different dual-tree wavelets 



Type of Dual- Tree Wavelets 


Q 


K 


Shannon wavelets (ideal) 


1 





B-spline wavelets, degree=l 


1 


0.9245 


B-spline wavelets, degree=3 


1 


0.0882 


B-spline wavelets, degree=6 


1 


0.0373 


Orthonormal spline wavelets, degree=l 


1 


0.9292 


Orthonormal spline wavelets, degree=3 


1 


0.1570 


Orthonormal spline wavelets, degree=6 


1 


0.0612 


Kingsbury's wavelets (q-shift Le Gall 5/3) [12] 


0.9992 


0.6586 



The practical challenge is the construction of different flavors of wavelets 
that fulfill or, at least, provide close approximations of these criteria. Indeed, 
the first criteria has been marked by an extensive research into the problem of 
designing both approximate and exact HT wavelet-pairs [16, 7]. 

We propose new design metrics for assessing the quality of the approxima- 
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tion. A simple measure for criterion (CI) is the correlation 



Q = 



max 



ll^^l||-||V'2|| 



, 



(20) 



The Cauchy-Schwarz inequality asserts that ^ g ^ 1 , where ^7 = 1 if and only 
if 'ip2{x) = Jif'ipi{x). Thus, the higher the value of p, the better would be the 
approximation. 

Next, note that criterion (C2) also has a simple Fourier domain characteri- 
zation: 



where g{u)) is complex- valued in general, and has a local axis of symmetry within 
its support. In particular, if g{uj) is constrained to be real-valued (correspond- 
ing to a symmetric ^){x)), then ^'1(2;) + j^2{x) has a constant instantaneous 
frequency (derivative of the phase) over its support. A reasonable quality metric 
for (C2) would then be the variation of the instantaneous frequency. Alterna- 
tively, we can also assess the degree of symmetry of S{co). In particular, we 
propose the following measure: 



That is, K, which lies between and 1, measures the disparity between S{(jj) 
and its reflection around the centroid. Indeed, k equal zero if and only if S{io) 
is symmetric (with Co as the centre of symmetry) . Conversely, a high value of k 
signifies greater local asymmetry in S{uj), and hence a poorer approximation of 
the modulation criterion. 

We computed the metrics g and k for different dual-tree wavelets (cf. Table 
1). The wavelets ipii^) and tp2{x) were synthesized using the iterated filter- 
bank algorithm, and the integrals involved in (20) and (21) were realized using 
high-precision numerical integration. It can easily be verified that g and k are 
appropriately normalized in the sense that they are invariant to the scale and 
translation of the synthesized wavelets. This is a necessary criteria since the 
wavelets are essentially synthesized by the filterbank algorithm at some arbi- 
trary scale and translation. 

The spline wavelets are analj^tic by construction [7], and hence ^) = 1 ir- 
respective of their degree. This is a necessary criteria since the synthesised 
wavelets are essentially at some arbitrary scale and translation. However, as 
their degree increases, the B-spline (resp. orthonormal) wavelets converge to 
the real and imaginary components of the complex Gabor (resp. Shannon) 




K = 



\S*{u) + uj)- S{ij - uj)\ dto 



(21) 



where u, the centroid of |5'(w)|, is specifled as 

_ ^ jr^to \S{u;)\ckv 
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wavelet which exhibits a symmetric spectrum (see §5 for details). The rapid 
decrease in the value of n reflects this improvement in symmetry (and also the 
rate of convergence). 

4 MULTI-DIMENSIONAL EXTENSION 

In this section, we extend the amplitude-phase representation derived in §3.1 to 
the multi-dimensional setting. The key ideas carry over directly, and the final 
expressions (cf. (27) and (16)) are as simple as their ID counterparts. The 
attractive feature of the multi-dimensional dual-tree wavelets is that, besides 
improving on the shift-invariance of the corresponding transform, they exhibit 
better directional selectivity than the conventional tensor-product (separable) 
wavelets [16]. For the sake of simplicity, and without the loss of generality, we 
derive the amplitude-phase representation for the particular two-dimensional 
setting. 

Two-Dimensional Dual- Tree Wavelets: To set up the wavelet notations, 
we briefly recall the construction framework proposed in [7] involving the tensor- 
products of one-dimensional analytic wavelets. Specifically, let <p{x) and (p'(a;) 
denote the scaling functions associated with the analytic wavelet t/ja{x) = ip{x) + 
jtp^x), where tp'{x) = Jfftp{x). The dual-tree construction then hinges on the 
identification of four separable multiresolutions of L^(R^) that are naturally 
associated with the two scaling functions: the approximation subspaces V{(p) (S 
V{(p),V{(p) (g) y((p'),y((p') (g) V{(p) and V{(p') (g) V{(p'), and their multiscale 
counterparts^. The corresponding separable wavelets - the 'low-high', 'high- 
low' and 'high- high' wavelets - are specified by 





= <p{x)i^{y), 


i)i{x) 


= (p{x)^'{y), 




= '>P{x)<p{y), 


■4>5{x) 


= i^{x)<p'{y), 


V'3(a;) 


= Hx)ij{y), 


•fpeix) 


= ^{xWiy), 




= v'{x)ij{y), 


i'loix) 


= V'ix)i;'iy), 


■tjjsix) 


= i^'{x)(p{y), 


■tjjnix) 


= ^'{xWiy), 


■tp^ix) 


= i>'{x)il){y), 


'4>i2{x) 


= ^'{x)^'{y), 



whereas the dual wavelets ■ipi,. . . ,'ipi2 are specified in terms of ip{x) and tp'{x) 
(here x = {x,y) denotes the planar coordinates). As far as the identification of 
the complex wavelets is concerned, the main issue is the poor directional selec- 
tivity of the 'high-high' wavelets along the diagonal directions. This problem 
can, however, be mitigated by appropriately exploiting the one-sided spectrum 
of the analytic wavelet tpaix), and, in effect, by appropriately combining the 

"^Thc tciisor-product V((p)^V{(p) denotes the subspELce spanned by the translated functions 
cp(- — m)cp(- — n), (m, n) S Z^. 
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wavelets in (22). In particular, the complex wavelets 

^'2(0;) = 'ilJa{x)(p'{y) = -ipsix) + j-ipii{x), 

"ifsix) = (p(x)V'a(y) = i'lix) + j-ipiix), 
*4(a;) = (p\x)ljja{y) = '4'7{x) + j4>lo{x), 

MX) = ^^rAx)My) = (^ ^^^-)^J^^^-) ) ^ ^a(.) -_^.(.) j (23) 

exhibit the desired directional selectivity along the primal orientations Oi = 
6I2 = 0, 6*3 = 6*4 = 7r/2, 6*5 = 7r/4, and 6»6 = 37r/4 respectively [7]. The dual 
complex wavelets ^i,...,^6) specified in an identical fashion using the dual 
wavelets i}p{x), are also oriented along the same set of directions. 

Directional Hilbert Transform (dHT) : Having recalled the complex wavelet 

definitions, we next recall the "quadrature" correspondence between the compo- 
nents of the complex wavelets that provides further insight into their directional 
selectivity. Akin to the HT correspondence, the components can be related 
through the directional HT ^,0 < 6* < tt: 

.3%f{x) ^ ~jsign{uju:)f{iv) 

where U0 = (cos 0, sin 0) denotes the unit vector^ along the direction 9. In 
particular, one has the correspondences 

Jm{^ e) =ng,M^f,) (£=!,..., 6), 

so that, by denoting the real component of the complex wavelet ^'^(a;) by tpe{x), 
we have the following convenient representations 

^e{x) = Mx)+jneMx) {£=!,..., 6) (24) 
which axe reminiscent of the ID analj^tic representation. 

4.1 Amplitude-Phcise Representation 

Let us denote the dilated-translated copies of the each of the six analysis wavelets 
^e{x) by §^,j,fe(a;), so that 

^e^i^kix) = Ei^k^e{x) (i G Z, fc € Z^), 

^Note that the half-spaces {a; : uJlj > 0} and {tti : uJlj < 0} play a role, analogous to 
that played by the half-lines ^ w ^ oo and — oo ^ w ^ in case of the HT, in specifying 
the action of dHT. 



15 



(f(x),i|Ji,k(x;0)) =21.1 



(a) 




(f(x),il)i,k(x;Ti[k])) =29.3 



(b) 

Figure 4: Wavelets corresponding to the step unit f{x): (a) reference wavelet 
'>Pi,k{x) = il'i,k{x; 0) corresponding to the conventional DWT; (b) shifted wavelet 
dual-tree V'i,fc(a;; The magnitudes of the signal correlation in either case 

clearly shows that the shifted wavelet has a better lock on the singularity. 
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where Si.fc is specified by 'E.i,k f{x) = 2'/(2*a; — fc). The corresponding dual-tree 
transform involves the analysis of a finite-energy signal f{x) in terms of the 
sequence of projections 

4[k] = \{f,^e,i,k), (25) 

where the use of the normalization factor 1/4 will be justified shortly. Before 

deriving the representation of f{x) in terms of the analysis coefficients c|[fe], we 
introduce the following fractional extensions of the dHT: 

= cos(7rr) y - sin(7rT) Jfe (r e R) (26) 

that formally allow us to capture the notion of direction-selective phase-shifts. 
Certain key properties of the fHT carry over to the fdHT. In particular, the 
family of fdHT operators {^/^.tItgr, ^ 6* < tt, exhibit the fundamental in- 
variances of 

• Translation: i^,r{/(- - y)}{x) = (i^,r/)(a; - y); 

• Dilation: 9{o,r{fiX-)}ix) = {9{e,rf)i\x); and 

• Norm: \\!Hs,rf\\ = \\f\\, for all / G L^{R^), 

which played a decisive role in establishing (11). Based on the above invariances, 
we can derive the following representation 

f{x)= J2 HH S.,fc{V^(x;r/[fc])} (27) 

involving the superposition of direction-selective synthesis wavelets affected with 

appropriate phase-shifts (cf. Appendix §7.3 for derivation details). In partic- 
ular, the wavelets ^e(x;T^[k]) are derived from the reference wavelet 'ipe{x) 
through the action of fdHT, corresponding to the direction 6e and shift r/[fe] = 
a,rg{cf[k]) /n As in the ID setting, further insight into the above representation 
is obtained by considering wavelets resembling windowed plane waves. 

4.2 Directional Modulated Wavelets 

A distinctive feature of the dHT is its phase-shift action in relation to plane- 
waves: it transforms the directional cosine cos(M^a;) into the directional sine 
sin{tigx). Moreover, what turns out to be even more relevant in the current 
context is that the above ac;tion is preserved for certain classes of windowed 
plane waves; in particular, we have 

Jfe {<p{x) cos{QuJx)} = (p{x) sm{QuJx) 

provided that ((>{x) bandlimited to the disk Dn = {w : ||a>|| < f2}. The following 

result - a specific multi-dimensional extension of Theorem 3.2 - then follows 
naturally for the fractional extensions: 
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Proposition 4.1 Let the window function (p(x) be bandlimited to the disk Dq. 

Then we have that 

iHo^r cos{flugx)} = (p{x) sm{Qu0X + ttt). (28) 

That is, the fdHT acts only on the phase of the oscillation while the window 
remains fixed. In particular, if the dual-tree wavelets are of the form 

ijjeix) = (pe{x)cos{CleU0^x) , £ = 1,...,6, (29) 

then the right-hand side of (27) assumes the form 

fixed window variable amp— phase directional wave 

53 Ve,i,k{^) ^hk[ \Ci[k]\cos{niuJ^x + TrTi[k]) |, 

{e,i,k) 

where (pe,i,k{^) the dilated-translated copies of the window (pe{x). The 
above expression explicitly highlights the role of r/ [fc] as a scale-dependent mea- 
sure of the local signal displacements along certain preferential directions. 

Indeed, this is the scenario for the spline-based transform proposed in [7] 
where the dual-tree wavelets asymptotically converge to the two-dimensional 
Gabor functions proposed by Daugman [8]. Moreover, these Gabor-like dual- 
tree wavelets were constructed using the B-spline scaling function and the semi- 
orthogonal B- spline wavelet. Replacing these with the orthonormal B-spline 
and the orthonormal wavelet, respectively, would then result in Shannon-like 
dual-tree wavelets - sinc-windowed directional plane waves - following the fact 
that the orthonormal spline multiresolution asymptotically converges to the 
Shannon multiresolution (cf. §3.2). flip the roles of the analysis and synthe- 
sis wavelets: we analyze the signal using the dual complex wavelet ^{x;a) = 
tp{x;a) + jil)'{x;a), and the Gabor-like wavelet ^'(a;;a) is used for reconstruc- 
tion. 



5 SHIFTABLE SPLINE WAVELETS 

If the wavelet is not modulated, we can characterize the action of the fHT and 
the fdHT by studying the family of wavelets {ip{x; r)},-eR and {ipe{x; r)}xeR) 
respectively. The remarkable fact is that it can be done explicitly for all spline 
wavelets derived from the fractional B-splines [5], which are the fractional ex- 
tensions of the polynomial B-splines. 

Spline Multiresolution: We recall that the B-spline P"{x), of degree a G 
M.Q and shift r e M, is specified via the Fourier transform 
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The degree primarily controls the width (and norm) of the function, whereas the 
shift influences the phase of the Fourier transform. As will be seen shortly, the 
latter property plays a key role in conjunction with the fHT. Though these func- 
tions are not compactly supported in general, their 0(l/|.x|"'^^) decay ensures 
their inclusion in L^(M) fl L^(R). More crucially, they satisfy certain technical 
criteria [5, 22] needed to generate a valid multiresolution of L^(R): 
(MRAl) Riesz-basis property: the subspace 

span^2{/?"(- - k)}kei. 

admits a stable Riesz basis. 

(MRA2) Two-scale relation; the refinement filter is specified by the transfer 
function 

(MRA3) Partition-of-unity property. 

Spline Wavelets: As far as the wavelet specification is concerned, the 
transfer function of the wavelet filter that generates the generic spline wavelet 
tp"{x), of degree a and shift r, is given by 

G«(e''^) = e'^Q°'{-e-^'')H^{-e-^", (31) 

where the filter Q°'{e^'^) satisfies the lowpass constraint Q"(l) = 1, and is inde- 
pendent of T. The filter plays a crucial role in differentiating between various or- 
thogonal (e.g., Battle-Lemarie wavelet) and biorthogonal (e.g., semi-orthogonal 
B-spline wavelet) flavors of spline wavelets of the same order [20, 21]. The as- 
sociated dual multiresolution is specifled by a dual spline function and a dual 
spline wavelet ■0"(a;); the fundamental requirement for the dual system is that 
the dual wavelet - generated using a transfer function similar to that in (31) 
- satisfies the biorthogonality criterion (V'"(' — "T')) V'r (' ~ ^)) = '^["^ ~ 
shall henceforth use the notation ip'^{x) to denote a spline wavelet of order a 
and shift r, irrespective of its genus (orthonormal, B-spline, dual spline, etc.). 

It turns out that the family of spline wavelets {V'" IreR (of a specified genus 
and order) is closed with respect to the fHT operation. 

Proposition 5.1 The fHT of a spline wavelet is a spline wavelet of same genus 
and order, but with a different shift. In particular, 

i^V"(a;) = C-r(a;). 

The above result (cf. part I of Appendix §7.4) signifies that the fHT acts only 
on the shift parameter of the spline wavelet while preserving its genus and order. 
Thus, for the dual-tree transform involving the corresponding wavelet basis 
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{V'ffe(2;)}(i,/s)ez2 and its HT pair, we have the following signal representation 

= ^ E,,k{\ci[k]\ 

involving the weighted sum of the appropriately "shifted" spline wavelets. 

Finally, we investigate the action of the fdHT on the 2D dual-tree wavelets 
ipe{x) constructed using a spline wavelet ■(/'"(a;) of a specific genus, and its 
HT pair V'r+i/2(^) I^l- turns out that, as in the ID case, the action is 
purely determined by the perturbation of the shift parameter of the constituent 
spline functions. However, the key difference is that the fdHT operators act 
"differentially" on the shifts of the spline functions along each dimension. Before 
stating the result we briefly digress to introduce a convenient notation. Observe 
that the six dual-tree wavelets are of the general form 

X] 9e{x; a, Tx)hi{y\ a, Ty) 

I 

where gi{x\ a, Tx) and ht{y\ a, Ty) arc spline scaling functions/ wavelets that have 
a common degree a but whose shifts depend on and Tj^ respectively. To 
explicitly emphasize the dependence on the parameters and Ty, we denote 
the dual-tree wavelets by ijji{x: a, t) with the shift- vector t = (t^, Ty) specifying 
the shift parameters of the spline functions involved along each dimension. For 
instance, the wavelets ipi{x;a,T) and tjj5{x;a,T) are specified (see (23)) by 

Mx;a,T)=i;l{x)p:^^{y), 

U^;a,T) = -^(C.(a;)C,(2/) -C.+i/2(^)C,+i/2(2/)), 

where t = (r, t) by construction. In general, setting r = (r, r) for all the six 
wavelets we deduce (see part II of Appendix §7.4 for a proof) the following: 

Proposition 5.2 The fdHT of a 2D dual-tree spline wavelet is a dual-tree spline 
wavelet of the same order and direction, but with a different shift: 

^e,ffpe{x; a,T) = (a;; a, r - f/x^tiej 

where yu^ = 1 for 1= 1, . . . , 4, and equals l/\/2 for i = 5 and 6. 

The result is quite intuitive. The horizontal and vertical wavelets can be 
'shifted' along the direction of the corresponding fdHT by perturbing the shift 
of the spline functions running along the same direction; the shift of the spline 
functions along the orthogonal direction remains unaffected. However, the di- 
agonal wavelets can be 'shifted' only by simultaneously by perturbing the shift 
of the splines along both dimensions. 
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Thus, as a direct consequence of (27) and proposition 5.2, we have the 
following signal representation: 

e,i,k 

where the shift information T^[k] = arg(cf[fc]), at different scales along each 
of the six directions, is directly encoded into the shift-parameter of the spline 
wavelets. As discussed in §4, for sufficiently large a, a, r) constructed us- 

ing the B-spline (orthonormal spline) wavelets resemble the Gabor (resp. Shan- 
non) wavelet where the shift T^[k] gets directly incorporated into the phase of 
the modulating plane wave. 

6 CONCLUDING REMARKS 

We derived an insight into the improved shift-invariance of the dual-tree complex 

wavelet transform based on single fundamental attribute of the same: the HT 
correspondence between the wavelet bases. Indeed, the identification of the 
fHT-transformed wavelets involved in representation (16) followed as a direct 
consequence of this correspondence. The shiftability of the transform was then 

established based on two key results: 

• the intrinsic invariances of the fHT group with respect to translations, 
dilations and norm-evaluations; and 

• theorem (3.2) describing the phase-shift action of the fHT on modulated 
wavelets. 

In particular, a multiscale ampUtude-phase signal representation was derived for 
the class of the modulated wavelets which highlighted the additional freedom 
of the wavelets to lock on to singularities of the signal. We also proposed 
certain metrics for accessing the modulation criteria and the quality of the HT 
correspondence between the dual-tree wavelets. These could prove useful in the 
design of new dual-tree wavelets with better shiftability. 

Before concluding, we would like to remark that the (direction-selective) 
shiftability of the dual-tree transform can also be extended to higher dimensions. 
In particular, the wavelet construction (23), the fHT correspondences (24), the 
modulation law (28), and, crucially, the amplitude-phase representation (27) 
carry over directly to the multi-dimensional setting. 

7 Appendix 

7.1 Proof of Theorem 3.1 

The sufficiency part of the theorem follows from the properties of the fHT 
operator; we only need to prove the converse. It is well-known that a unitary 
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linear operator T on L^(R) is translation invariant if and only if there exists a 
bounded (complex-valued) function m{uj) such that 

ff{uj)=m{oj)f{u), 

for all / G L^(R.) [18, Chapter 1]. This Fourier domain characterization reduces 
the problem to one of specifying a bounded function m(a;) such that T has the 
desired invariances. It can be readily demonstrated that the dilation-invariance 
criterion translates into the constraint 

m{au) = m{u}), a > 0. (32) 

Moreover, since the real and imaginary components of to(w) must independently 
satisfy (32), one even has TOi(aw) = mi(cj) and m2(aa;) = m2(w), for all a > 0, 
where mi{uj) and m2(w) are the real and imaginary components of m{oj). 
Next, observe that the Hermitian symmetry requirement 

m*{uj) = m{-(j) (33) 

on the multiplier require mi(w) and m2(w) to be even and odd symmetricre- 
spectively; that is, mi(— w) = mi{ui) and m2{—cjj) = — m2(cj). These constitute 
the crucial relations, since one can easily verify that the only bounded functions 
(up to a scalar multiple) that satisfy (32) and (33) simultaneously are the con- 
stant function mi(u;) = 1, and its "skewed" counterpart 7712(0;) = sign(a;); that 
is, it is both necessary and sufficient that 

m{u) = 71 + i72 sign(a;), 

for some real 71 and 72. Finally, combining the equivalence 

l|T/|P = ^ / M'^)n/HPdc.= (7i'+72')'ll/ll' 

obtained through Parse val's identity, with the norm invariance requirement we 
arrive at the criterion 7i -I-7I = 1. Therefore, it is both necessary and sufficient 
that m{u) = cosO + j sm9 sign(a;) for some 6 G (—• tt, tt], and this establishes 
the representation 

T = cose y -sine Jif. 

The uniqueness of ^ e (— 7r,7r] necessarily follows from the quadrature corre- 
spondence in (6). □ 

7.2 Proof of Theorem 3.2 

We note that the Fourier transform of f{x)g{x) is given by the convolution 
{2'k)~^ {f * g){uj) . Thus, if we denote the Fourier transform of Hr\^f{x)g{x)^ by 
F{(jj), then following definition (7), we have that 
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In particular, this gives us the Fourier representation 

^{f{x)g{x)] ^ ^ J^F{co)e^--du; 

^ - duj. (34) 

Now, by commuting the order of the integrals and by applying the frequency 
translation ( = co — ^, we can rewrite (34) as the double integral 

K I eJ'^^/(ae^''^'^'8"('^) 5(w-0dwd^ 

^ / e^«+^)^/(0 e^'^^^'s^Ci+O giQdCdi. (35) 

Jo 
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The effective domain of integration in (35) gets restricted to the region Q = 
{{Ci ■ ICI ^ 1^1 < ri} as a consequence of the assumptions on the supports 
of /(^) and g{C)- In particular, it can cfisily be verified that sign(^ + i^) = sign(^) 
on Q; this allows us to factor (35) into two separate integrals giving the desired 
result: 



= f{x)Hrg{x). □ 
7.3 Derivation of Equation (27) 

As a first step, we consider the equivalent representations of f{x) in terms of the 
four distinct bases generated through the dilations-translations of the separable 
wavelets in (22): 

(i,k) 

(36) 

for p = 0, 3, 6 and 9, where the expansion coefficients are specified by 

= (9=1,..., 12). (37) 



23 



Next, we combine and regroup these expansions as follows 

fi^) = m W ("2,i[fe]^2,i,fc(a;) + a8,i[fe]V'8,i,fc(a;)) + ^(a5^i[k]'ip5^i^k{x) + aii^ilkj-ipii^i^kix) 
(j,fe) I 

+ ^(«i,i[^]V'i,j,[fc](a;) + a4,4fc]'04,i,fe(ic)) + |(a7,i[fc]V'7,i,fc(a;) + aio,j[fc]V'io,i,fc(a;)) 

, 1 r ri,i r.n /V'3,»,fc(a;) - i/'i2,j,fc(a5)\ , 1 . r. /■i/'6,i,fc(«) + V'9 
4:71 - ai2,i[k]) l-^ ^ 1 + 4:^('*6,4fc] + ag.iW) {— ^ — 

, 1 ^ ri,i , r,n /^V'3,i,fe(a;) + V'i2,i,fe(a;)\ 1 , , , / i^Q^i^kix) - tpg 



The terms on the right-hand side can now be conveniently expressed in terms 

of the cocfRcicnts cf[fe] = |cf [fe] |c^'^»['^l and the dual-tree wavelets i^e.i,k{x). For 
instance, consider the terms in the third line. The wavelet pair ip5,i,k{x) and 
^54'5,i,k{x) are readily identified; moreover, the correspondences 

= 'fMf,^..,k) = -\c![k]\sinm, 

that follow from (25) and (37), allow us to rewrite it as 

|cf[fc]|(cos<^f[fc]V5,i,fc(a?) -sin<^f[fc]:?/e5^5,i,fc(a;)) = |cf[fc]| ^p5,i,k{x;4[k]) 

with the shift specified by T^[k] = [fe]/7r. Subjecting the rest of the terms to 
a similar treatment we arrive at the desired representation. 



7.4 Proof of Propositions 5.1 and 5.2 

Part I: This result is easily established using an auxiliary operator. Specifically, 
we consider the fractional finite-difference (FD) operator 

Arf{x)^Dr{e^nf{Lj), (38) 

corresponding to the frequency response D"(e-''^) = (l — e^^")^ (l — e^'^) ^, 
which allows us to relate fractional B-splines (and the corresponding filters) of 
the same order but with different shifts. In particular, it allows us to express 
the action of the fHT on a B-spline as a linear (digital) filtering: 

5^/?^(x) = A^/3^_^(x) 

= Y,Mk]Pr-f{^-k). (39) 
feez 
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Indeed, based on definitions (30) and (38), and the identity^ (juj) '^(—juY = 
exp ( — j7rTsign(a;)), we have the following factorization 



/3^?(a;) = (ia;)-^(-ia;)^D^(e'--)/3«_^(w) 

= exp ( - jTTf sign(a;))DT=(e-"^)/3"_^(a;) 

which results in the equivalence 

:Hfl3'^{x) ^ exp (jttt sign(c^))/3;?(a;) 
= D,(e^'-)/3."-,(^) 
^ A,/?°_,(x). 

Next, we observe that the conjugate-mirrrored version of the FD filter can 
also be used to relate the spline refinement (scaling) filters of the same order but 

with different shifts: H^_^{e^'^) = Dfi-e'^'^) i?^(e^"). In particular, wo have 
the relation g"_f[k] = {df * 5") [A:] between the corresponding wavelet filters: 

G«_^(e-'''^) = e^''Q'^{-e-^'')H^_^{-e-^^) 

= D^(e^")e^"Q"(-e-^'^)if^(-e-^''") 

= D^(e^-)G^_^(e^-). (40) 

As a result of (39) and (40), we have the desired equivalence: 

fcez 

= J29r[k]{ E Mn]P?-A- - n)}{x - k) 

fcez nez 

= ^idr * df)[m]Pr-f{x - m) 
mez 

= ^lj^_,{x/2). 

Part II: We derive the relation for the spline wavelets ipi{x;a,T) and 
ip5{x; a,T) (the rest can be derived identically). Using proposition 5.1, we im- 
mediately arrive at one of the results: iN^^^f'4'i{x;a,T) — .Wb,f {V'" (•^)}/^t (y) = 

The second result relies on the factorization i?4/4,T = ^,t/2^/2,t/2 that 
holds for functions whose frequency supports are restricted to the quadrants 
{{ujx,'-0y) : tOx > OjUjy > 0} and {{uj^jUty) : w^. < OjWj, < 0}. In particular, the 
condition is satisfied by ip^iu); a, r) so that, in conjunction with proposition 5.1, 

^We specify the fractional power of a complex number 2 by z'' = |2|TeJ'i"i''s(2) correspond- 
ing to the principal argument | arg(2;)| < n. On this branch, the identity (2122)'*' = zjz2 holds 
only if aj:g(2i) + arg(2i) G (— 7r,7r) [17, Chapter 3]. 
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we have 
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